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摘 要 : 叶绿素 含量 是 表征 植被 生长 状况 的 重要 参考 指标 ,利用 高 光谱 技术 快速 ,精确 地 监测 棉花 叶片 叶绿素 含 


量 , 以 新 疆 125 个 苗 其 


i 


棉花 叶片 样本 为 研究 对 象 ,通过 测定 其 叶绿素 含量 与 光谱 数据 ,采用 多 种 光谱 预 处 理 和 多 植 


被 指数 相 结 合 的 方法 ,构建 了 WOA-RFR 棉花 叶片 叶绿素 含量 定量 反 演 模型 ,并 与 SVR 和 RFR 模 型 结果 进行 对 比 
分 析 。 结 果 表 明 :(1) 光谱 变换 方法 中 对 数 变 换 分数 阶 微分 和 连续 小 波 变 换 均 能 有 效 地 提高 植被 指数 与 叶绿素 含 


量 的 相关 性 。(2) 基于 分 数 阶 微分 0.9 阶 变 


a 


的 Vogelmann3 , RVI, DVI SR {67s-100) 7, Mndvizs, ND , VOG1 , NVI, TVI $H 


V0G2 植 被 指数 组 合 的 WOA-RFR 模型 反 演 效果 最 佳 ,其 建 模 集 和 验证 集 模型 尼 分 别 为 0.920 和 0.95$5,RMSE 分 别 


H 0.987 710.986 , MRE 分 别 为 0.013 和 0.014 ,与 RFR 和 SVR 模型 相 比 ,预测 精度 有 所 提高 ,WOA 算 法 优化 模型 效果 


明显 。 研 究 结果 可 为 棉花 叶片 叶绿素 含量 定量 反 演 提供 决策 依据 。 
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叶绿素 含量 是 表征 作物 健康 状况 的 重要 指标 ， 
其 与 作物 生长 紧密 相关 "。 叶 绿 素 含量 的 传统 测定 
方法 成 本 高 ,操作 复杂 效率 低 , 难 以 满足 精准 农业 
快速 发 展 的 需求 ”。 目 前 ,精确 .快速 ,高效 的 高 光 
谱 技 术 被 证 明 在 植被 叶绿素 含量 获取 方面 有 一 定 
的 优势 ”。 当 前 棉花 叶片 叶绿素 含量 的 光谱 反 演 主 
要 采用 多 种 预 处 理 方式 相 结 合 的 方法 ,分析 光谱 与 
叶绿素 含量 的 相关 性 ,并 选择 特征 波长 或 者 构建 高 
光谱 植被 指数 ,建立 叶绿素 含量 监测 模型 “”。 前 
人 利用 多 种 光谱 预 处 理 方法 ,如 小 波 变 换 汪 一 阶 导 
数 卸 二 阶 导数 “分数 阶 微分 变换 "等 ,有 效 减 少 其 
至 消除 高 光谱 的 元 余 信息 ,选择 合适 的 光谱 变换 方 
法 能 提升 作物 叶绿素 反 演 精度 。 同 样 ,植被 指数 法 
作为 高 光谱 反 演 农作物 生理 生化 指标 的 重要 方法 ， 
也 已 被 广泛 应 用 于 反 演 作物 的 生物 化 学 参量 。 
随 着 机 器 学 习 算 法 快速 发 展 , 多 元 逐步 回归 (Multi- 
ple Linear Regression, MLR) 2 、 偏 最 小 二 乘 回 归 
(Partial Least Squares Regression , PLSR )52 随机 森林 
回归 (Random Forests Regression, RFR)‘. 支持 向 量 
机 回归 (Support Vector Machine Regression, SVR) | 
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叶绿素 含量 ;鲸鱼 优化 算法 


等 方法 在 估算 叶片 叶绿素 含量 研究 中 均 取 得 了 较 
好 的 效果 。 茹 克 亚 等 “利用 对 数 变换 和 分 数 阶 1.2 
阶 微分 计算 的 植被 指数 (NDVI、DVI、RVI、MSR;、 
MSRooso、C1) 建立 的 PLSR 模型 的 估算 精度 最 优 。 
Liu 等 利用 连续 小 波 变换 对 马 铃 暮 叶 绿 素 含 量 进 
行 分 析 , 结 果 表 明 连 续 小 波 变换 方法 能 够 有 效 的 提 
取出 对 马铃薯 叶绿素 含量 较为 敏感 且 稳 健 的 小 波 
特征 ,增强 了 光谱 与 叶绿素 含量 之 间 的 相关 性 。 张 
思 楠 等 "通过 植被 指数 对 多 枝 梭 柳 叶片 叶绿素 含量 
估算 研究 , 选 出 了 结构 简单 且 表 现 较 好 的 光谱 指 
数 ,以 上 研究 表明 ,植被 指数 法 在 反 演 作物 叶片 生 
理 生 化 指标 方面 ,能 够 取得 较 好 的 效果 。 而 郭 超 凡 
等 中 提出 的 多 植被 指数 协同 估算 模型 精度 相对 于 
单 植 被 指数 模型 精度 明显 提高 。Qi 等 对 花生 叶 
片 叶 绿 素 含量 的 高 光谱 反 演 结果 表明 ,NDSI、RSI、 
DSI、SAVI 是 最 佳 的 光谱 指数 。 众 多 学 者 利用 不 同 
机 器 学 习 算 法 ,不 同 光谱 预 处 理 方法 ,以 及 植被 指 
数 法 对 作物 属性 的 定量 监测 进行 了 大 量 的 研究 并 
取得 了 显著 的 成 果 。 因 此 ,选择 合适 的 光谱 预 处 理 
方法 .植被 指数 法 与 机 器 学 习 方 法 ,对 于 利用 光谱 
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反 演 作物 叶片 叶绿素 含量 尤为 重要 。 

综 上 所 述 ,光谱 数据 变换 方法 和 植被 指数 法 广 
泛 应 用 于 作物 参量 佑 测 研究 中 。 但 多 数 研 究 均 是 
基于 单一 的 光谱 数据 变换 或 单一 植被 指数 构建 模 
型 ,而 多 种 光谱 数据 变换 和 多 植被 指数 组 合 佑 测 棉 
花 叶 绿 素 含量 的 研究 鲜 有 报道 。 因 此 ,本文 结 合 前 


南 划 ,塔里木 盆地 中 北部 ,是 典型 的 而 形 平 原 绿 
洲 。 地 形 北 高 南 低 , 自 西北 向 东南 倾斜 ,最 高 海拔 
高 程 为 4550 m ,最 低 海 拔高 程 922 m。 地 处 暖 温带 ， 
气候 干燥 ,降水 稀少 ,夏季 炎热 ,冬季 干 冷 ,年 温差 和 
日 温差 均 较 大 , 属 暖 温带 大 陆 性 干旱 气候 ,年 均 温 为 
10.5~11.4 ,年 均 降 水 量 范围 为 43.1~51.6 mm", JE 


人 研究 选择 46 种 在 农作物 参量 研究 方面 应 用 较 广 
的 植被 指数 ,基于 传统 数学 变换 ,分数 阶 变换 ,小 波 
变换 ,结合 相关 系数 法 与 SVR、RFR 和 鲸鱼 优化 算 
法 (Whale Optimization Algorithm, WOA ) 优 化 RFR 建 
立 多 种 植被 指数 回归 模型 , 旨 在 解决 以 下 问题 : (1) 
光谱 数据 变换 方法 是 否 能 提高 植被 指数 与 棉花 叶 
片 叶绿素 含量 的 相关 性 。(2) 探究 最 佳 的 植被 指数 
组 合 和 最 佳 估 测 模型 。((3) WOA-RFR 估 测 模型 较 
RFR 模型 的 优化 效果 。 


1 材料 与 方法 


1.1 研究 区 概况 
人 研究 区 位 于 新 疆 阿 克 苏 地 区 库 车 市 阿 克 吾 斯 
HË S (82°58 ~83°30'E, 41°16’~41°31'N) ,地 处 天 山 


(a) 渭 库 绿洲 
82°30'0"E 


83°0'0"E 


(b) Bal ot 


83°0'0"E 


41°30'0"N 


41°0'0"N 


83°0'0"E 83°30'0"E 


(©) 棉花 冠 层 


(d) 叶绿素 计 


于 干 星 与 极端 干旱 区 。 研 究 区 及 样 点 分 布 情况 见 
图 1。 
1.2 棉花 样本 采集 和 叶绿素 含量 的 测定 

棉花 样本 数据 的 采集 于 2021 年 5 月 26 日 至 6 
月 3 日 ,在 阿克苏 地 区 库 车 市 阿 克 吾 斯 塘 乡 的 农用 
棉田 进行 ,棉田 总 面积 为 0.5 hm2, 棉 花 品 种 是 新 陆 
中 35 号 ,棉花 处 于 苗 期 阶段 , 共 采 集 125 个 样本 ,采样 
点 分 布 如 图 1 所 示 , 按 照 行距 为 70 cm, 株距 为 10 cm 
的 地 膜 覆 盖 礁 种 进行 种 植 ,棉田 的 施肥 ,栽培 管理 ， 
按照 当地 种 植 要 求 正常 进行 。 

选用 叶绿素 计 (SPAD-502Plus ) 在 获取 高 光谱 
数据 的 同时 测定 苗 期 析 花 叶片 SPAD 值 ,叶绿素 计 
(SPAD-502Plus ) 测 得 的 读数 与 叶绿素 含量 有 密切 联 
系 , 可 以 将 SPAD 值 转换 为 叶绿素 含量 进行 估 测 ” ， 
因此 ,本 文 将 SPAD 值 等 同 于 叶绿素 含量 进行 研 


HES 
(c) 棉花 采样 点 
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图 1 研究 区 及 样 点 分 布 示 意图 


Fig. 1 Schematic diagram of the distribution of the study area and sample points 
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究 。 测 量 每 个 样 点 上 中、 下 3 层 共 6 片 叶 子 ,测量 
每 片 叶子 的 不 同 部 位 6 次 , 取 其 平均 值 作为 最 终 值 。 
1.3 光谱 数据 采集 与 处 理 

棉花 冠 层 光谱 反射 率 使 用 ASD Field Spec 
Hand Held 便携 式 光 谱 仪 测定 ,该 仪器 的 探测 波长 
范围 为 325~1075 nm, 光谱 分 辩 率 为 3 nm。 光 谱 测 
量 时 间 为 北京 时 间 11:30 一 15:30, 在 无 风 无 云 晴 朗 
的 天 气 下 进行 。 测 量 前 先进 行 白 板 校 准 , 测 量 时 棉 
布 冠 层 与 光谱 仪 探 尖 的 距离 为 25 cm, Ff HE ELT] 
下 ,光谱 仪 的 视 场 角 为 25° ,光谱 扫描 时 间 至 少 为 8 s， 
每 个 采样 点 重复 测量 6 次 ,每 测量 3 次 后 进行 一 次 
白板 校准 。 最 后 将 在 ViewSpec Pro 软件 中 计算 得 出 
的 每 个 采样 点 的 6 条 曲线 的 平均 值 作为 各 采样 点 的 
光谱 反射 值 。 

本 文 对 原 波 段 光 谱 反 射 率 (Original Reflec- 
tance, OR) 进行 对 数 变 换 (Logarithmic Transforma- 
tion Spectrum, LogR ) 、 平 方 根 变换 (Square Root 
Transformation, JR ) 、 包 络 线 去 除 (Continuum-re- 
moval Transformation, CR) 、 分 数 阶 微分 变换 (Frac- 
tional-order Difference, FD ) .连续 小 波 变换 (Continu- 
ous Wavelet Transform, CWT) 等 ,其 中 传统 的 对 数 变 
换 ,平方根 变 换 能 够 使 可 见 光 区 域 的 光谱 差异 增 
强 , 使 光谱 的 吸收 ` 反射 特征 更 加 突出 ; 包 络 线 去 除 
处 理 可 压抑 背景 光谱 ,并 扩大 弱 吸 收 特征 信息 "; 分 
数 阶 微分 变换 有 增多 敏感 波段 以 及 提高 相关 系数 
的 作用 “ ;连续 小 波 变 换 能 够 挖掘 被 隐藏 有 效 光 谱 ， 
提升 光谱 信息 ,增强 光谱 特征 。 

1.4 植被 指数 的 提取 

本 文 结合 前 人 的 研究 “” ,选取 广泛 应 用 于 农 
作物 叶绿素 含量 研究 的 46 种 植被 指数 (VI) ,作为 叶 
绿 素 含 量 反 演 模型 的 参数 ,如 表 1 所 示 。 

1.5 模型 构建 和 评价 

本 研究 利用 多 元 回归 的 支持 向 量 机 回归 (SVR ) 
模型 和 随机 森林 回归 (RFR ) 模 型 分 别 进 行 建 模 一 ， 
比较 其 建 模 效 果 ,并且 采用 鲸鱼 算法 (Whale Optimi- 
zation Algorithm, WOA ) 优 化 RFR 模 型 ,检验 模型 优 
化 效果 ,最 终 往 选 出 建 模 效果 最 好 的 模型 。 将 125 
个 样本 分 为 两 组 ,一 组 为 建 模 集 样本 (n=82) 用 于 模 
型 的 建立 ,一 组 为 验证 集 样本 (n=43) 用 于 模型 的 
验证 。 

模型 的 精度 主要 由 决定 系数 (Determination Co- 
efficients, R’) 、 均 方 根 误差 (Root Mean Squared Er- 


ror, RMSE) 和 平均 相对 误差 (Mean Relative Error, 
MRE) =i DEVE ff, RMSE 和 MRE 越 小 ,说 明 
模型 的 误差 小 , 尼 越 接近 于 1 ,说 明 模型 的 拟 合 效果 
越 好 。 

鲸鱼 优化 算法 是 2016 年 提出 来 的 一 种 新 型 智 
能 优化 算法 ,其 源 自 于 对 自然 界 中 座 头 鲸 群体 狩猎 
行为 的 模拟 ,通过 鲸鱼 群体 搜索 包围. 追捕、 攻击 
猎物 等 过 程 实现 优化 搜索 的 目的 ”。 该 算法 的 计 


算 方 程 及 原理 如 下 : 
WOA 寻 优 路 径 更 新 方法 为 : 
D=ICX'(t) -XGO (1) 
X(t+1)=X'(t)-A-D (2) 


式 中 : X 表示 最 优 解 坐标 ; 工 为 当前 解 坐标 ,系数 
向 量 4 和 C 分 别 为 


A=2a:r-a (3) 
C=2r (4) 
式 中 :7 为 [0,1] 的 随机 向 量 。 
Be — JAB SK Sk MAS AT X ang ,搜索 猎物 数学 公式 
如 下 : 
D=IC-X,,,, -X (5) 
X(t+1)=X,,,.4-A°D (6) 
螺旋 攻击 方法 如 公式 47) 所 示 : 
XU+I=D'e“…cos(2m0)+XOO (7) 
当 座 头 鲸 发 现 猎物 后 ,进行 包围 捕食 或 螺旋 攻 
击 ,具体 由 选择 概率 决定 。 


X(t+1)= aS mas 
~ | Deel’ cos(271) + X"(t) 


p<0.5 
p=0.5 

(8) 
通过 不 断 更 新 YU+D ,按照 公式 (8) 进 行 猎 食 
运动 ,循环 迭代 直至 获取 最 优 个 体 。 


2 结果 与 分 析 


2.1 棉花 叶片 叶绿素 含量 统计 与 光谱 特征 分 析 
建 模 集 和 验证 集中 棉花 叶片 叶绿素 含量 的 平 
均值 分 别 为 56.74 和 54.98 ,标准 差分 别 为 2.49 和 
2.73 ,变异 系数 分 别 为 4.4% 和 5.0%; 全 部 样本 的 平均 
值 为 56.13 ,标准 差 为 2.70, 变 异 系数 为 4.8% ,三 个 数 
据 集 变异 系数 均 小 于 109% ,数据 变异 性 小 (图 2)。 

为 分 析 棉 花 叶片 叶绿素 含量 与 光谱 反射 率 的 
响应 特征 ,计算 全 部 样本 叶绿素 含量 的 较 高 值 BE 
低 值 .平均 值 和 中 间 值 ,研究 了 不 同 叶绿素 含量 水 
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表 1 高 光谱 植被 指数 
Tab. 1 Previously published hyperspectral vegetation indices 
植被 指数 计算 公式 文献 植被 指数 计算 公式 文献 
GNDVI (Roo = Rsoo/(Reoo + Rsoo) 19] NRI (Rsn = Re)/(Rsn + Rezo) 20] 
PSRI (Regs = Rsoo)/Ras0 20] NPCI (Reso = Rino) (Reso + Razo) 20] 
VARI (Rsss— Reso)/(Rsso + Reso + Raso) 21] RVI Ril Reo 21] 
NDVI Reo Reso! Rsoo + Reso 21] DVI Ra — Reco 21] 
voGl Ryo/Roao 21] VOG2 Roy Ray/Rns — Row 21] 
fel Rw/R 23] Lichtenthalerl 及 ao/Rem 23] 
Lichtenthaler2 (Roo — Reso)/(Reoo + Reso) 23] SIPI (Reco — Raso/(Reo + Riso) 23] 
PSSRa Roo! Reso 23] Vogelmann1 (Rna = Rr)/(Rns + Rro) 23] 
Vogelmann3 Rigo! 23] ND (Ross = Rios)/(Rsss + Roos) 23] 
Clred edge RaRa 24] Clgreen Roo! Rsso-1 24] 
RI-half Ry/Rios 24] GRVI (Rero/Rsoo) — 1 24] 
RDVI (Rsoo — Rezo)! {Reo ~ Ren 25] Dattl (Resp — Roig) (Reso + Reso) 25] 
Datt2 Rs /Ro 25] Datt3 R/R 25] 
Cartel Res /Rs 25] Carte2 Roos! Roy 25] 
Carte3 Roos! Rasa 25] Carte4 Rol Roe 25] 
SR 1325600) R/R 25] SRiszs70 RalR 25] 
NVI (Rj; — Ray)/Reos 25] MSR705 (Raiso ~ Rass)! (Ras + Raas) 25] 
SAVI (1 + 0.5) X(Rgog — Rero)/(Reoo + Ren + 0.5) 19] 
GARI [Reoo ~[Rsso — 1.7 X (Rago — Rezo) Reon +LRsso — 1-7 X (Ri + Revo) 19] 
OSAVI 1.16 X[(Rs — Reon)/(Rsoo + Ren +0.16) ] 21] 
Ro 红 谷 (640 ~ 680 nm) 内 反射 率 最 小 值 22] 
mNDos (R — Raps)!(Reso + Rags — 2 X Rys) 24] 
EVI 2.5 X[(Rgoy — Rero)/(Bgoy + Rony -7.5 X Rigs + 1) 25] 
MCARI [(Riw = Row) = 0.2 X (Rr = Rsso)] X (RoRo) 25] 
TVI 0.5 [120 X (Rs — Rsso) ~ 2-5 X (Reno — Rsso)] 25] 
MTVI1 1.2 X [1.2 X (Rego — Rsso) — 2-5 x (Re — Rsso) 25] 
REP 700 + 40 X [(Rer + Royo )/2 = Ri (Ri — Row) 25] 
SPVI 0.43.7 X (Roa — Reso) — 1-2 | Rz- Rool 25] 
SPVI2 0.4 X 3.7 X(Reoy — Rego) — 1-2 | Reso — Reso 25] 
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平 下 的 原始 光谱 反射 率 ( 图 3)。 不 同 叶绿素 含量 样 
本 的 光谱 曲线 的 变化 大 体 相似 但 也 有 一 定 的 差别 ， 
这 表明 叶绿素 含量 的 多 少 对 光谱 反射 率 有 一 定 的 影 
响 ,在 可 见 光波 段 范围 内 550 nm 附近 形成 了 一 个 反 
射 峰 ,690 nm 附近 形成 了 一 个 吸收 和 谷 ,700~750 nm 
波段 范围 内 ,反射 率 急 剧 上 升 直 至 形成 一 个 陡坡 ， 
950 nm 处 反射 率 逐 渐 下 降 ,至 975 nm 又 逐步 上 升 。 
325~680 nm 波段 范围 内 叶绿素 含量 低 反 射 率 越 高 ， 
这 也 证 明了 棉花 叶片 叶绿素 含量 对 可 见 光 波段 的 
敏感 性 。 

由 图 4 可 知 ,对 数 变 换 和 包 络 线 统 去 除 变 换 的 
相关 系数 在 全 波段 没有 提高 ,相关 性 低 于 原始 光 


谱 , 分 数 阶 变换 的 相关 系数 (以 FD-1.3 为 例 ) 在 700~ 
800 nm 相关 性 达到 最 高 (r=0.41) ,小 波 变 换 ( 以 
CWT-7 为 例 ) 的 相关 性 在 600~700 nm 呈 负 相关 
(r= -0.30) ,在 900~1000 nm 呈正 相关 ,相关 系数 最 
高 达到 (r=0.34) 。 对 比 不 同 数 学 光谱 变换 方法 的 
Pearson 相关 系数 ,分 数 阶 微分 变换 和 小 波 变换 是 较 
好 的 预 处 理 方法 。 
2.2 植被 指数 组 合 和 筛选 

本 文 运用 46 种 植被 指数 ,基于 原始 光谱 、 包 络 
RER .分 数 阶 变换 (以 0.3 为 步 长 ) .小 波 变换 (分 
EREL 23 2°27 2) 等 处 理 的 光谱 数据 上 进行 植 
被 指数 运算 , 选 出 若干 个 通过 0.01 相关 性 显著 水 平 
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注 :SPAD 为 相对 叶绿素 含量 ;SD 为 标准 差 ;C, 为 变异 系数 。 
图 2 棉花 叶片 叶绿素 含量 统计 


Fig. 2 Statistics on chlorophyll content in cotton leaves 
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图 3 不 同 叶 绿 素 含量 下 光谱 反射 率 变化 


Fig.3 Change in spectral reflectance under different 


chlorophyll content 


的 植被 指数 ,作为 一 个 组 合 植被 指数 (图 5)。 原 始 
光谱 数据 上 通过 0.01 相关 性 水 平 的 植被 指数 只 有 6 
个 ,而 传统 数学 变换 方法 中 的 对 数 变换 ,分数 阶 变 
换 中 的 0.3 .0.6.0.9 阶 ,小 波 变换 中 的 第 3、 第 7 尺度 
中 通过 0.01 相关 性 显著 水 平 的 植被 指数 数量 超过 
10 个 ,说 明 这 些 数 学 变换 形式 提高 了 更 多 植被 指数 
与 棉花 叶片 叶绿素 含量 的 相关 性 ,对 这 些 变换 上 的 
植被 指数 按照 与 叶绿素 含量 的 相关 性 系数 大 小 排 
序 , 取 相关 系数 大 的 前 10 个 植被 指数 作为 一 个 组 合 
植被 指数 ,并 且 选 择 这 些 植被 指数 组 合 进行 下 一 步 
分 析 。 在 平方 根 变换 、 包 络 线 变换 分数 阶 的 1.2 
阶 、 小 波 变换 的 第 1、 第 5 和 第 9 尺度 变换 形式 上 通 
过 0.01 相关 性 显著 水 平 的 植被 指数 个 数 较 少 ,分数 
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注 :FD-0.3 为 分 数 阶 1.3 阶 变换 ;OR 为 原始 光谱 反射 率 ;LOG 为 对 数 
光谱 变换 ;CR 为 包 络 线 去除 ;CWT-7 为 小 波 变 换 第 7 个 尺度 ,下 同 。 
图 4 不 同 数学 变换 反射 率 与 叶绿素 含量 的 相关 性 分 析 


Fig. 4 Correlation analysis of reflectance and chlorophyll 


content in different mathematical transformations 


阶 1.5 阶 及 之 后 的 阶 次 上 通过 0.01 相关 性 显著 水 平 
的 植被 指数 个 数 均 少 于 1 ,不 做 进一步 分 析 。 
2.3 多 植被 指数 组 合 相 关 性 分 析 

不 同 光 谱 变 换 下 植被 指数 与 棉花 叶绿素 含量 
相关 性 统计 表 , 如 表 2 所 示 ,每 一 种 光谱 变换 下 的 植 
被 指数 与 棉花 叶片 叶绿素 含量 之 间 进 行 相关 性 分 
析 , 原 始 光 谱 下 植被 指数 的 相关 性 均 低 于 光谱 变换 
后 植被 指数 的 相关 性 ,表明 对 光谱 数据 进行 不 同 数 
学 变换 ,能够 提高 植被 指数 与 叶绿素 含量 的 相关 
性 ,分 数 阶 0.6 阶 次 下 植被 指数 的 相关 性 为 -0.33~ 
0.34 ,小波 变 换 的 第 7 尺度 下 则 为 -0.31~0.34, 证明 
分 数 阶 变换 和 小 波 变换 更 能 提高 叶绿素 含量 和 植 
被 指数 相关 性 。 
2.4 建 模 结果 与 分 析 

利用 多 种 光谱 预 处 理 方 法 和 不 同 植被 指数 组 
合 ,分 别 建立 SVR RER 以 及 WOA-RFR 模型 估算 棉 
花 叶 片 叶绿素 含量 ,研究 不 同 植被 指数 组 合 方法 对 
模型 精度 的 影响 ,其 结果 如 表 3 所 示 。 从 表 中 可 知 ， 
在 不 同 数 学 变换 上 RFR WOA-RFR 模型 效果 较 好 ， 
且 WOA-RFR 优 于 RFR 模 型 ,证 明了 鲸鱼 优化 算法 
的 可 行 性 ,而 SVR 模 型 的 反 演 效果 相对 较 差 ,SVR 
模型 在 建 模 集 和 验证 集 的 模型 精度 差异 较 大 ,这 说 
明 该 模型 受 样本 个 数 影响 较 大 ,模型 稳定 性 差 。 基 
于 对 数 变 换 多 植被 指数 组 合 的 RFR 和 WOA-RFR 模 
型 的 训练 集 尼 分 别 为 0.884 和 0.912 ,验证 集 尼 分 别 
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10 f VOG1 SAVI DVI DVI Lichtenthaler2 REP 
NDCI OSARI SPVI2 ND CIred MCARI 
8r Datt2 DVI SAVI Vogelmann3 PSSRa Datt2 Vogelmannl 
x ND ND OSARI SRres70 Vogelmann3 MSR7s ”SRpsz-eol Dattl 
< 6F SAVI RVI RDVI EVI RVI RVI RVI SPVI GRVI 
<a DVI NDVI Carte4 ND mNDs MSRyos VOG1 DVI MTVIl 
pesn 4F ND GNDVI SPVI RDVI VOG1 mND ps5 SIPI PSRI cartel Lichtenthaler1 
RDVI PSSRa NPCI SPVI2 Carte? NVI SRyers.700) Clgreen ND Ro SAVI 
27 SPVI ClIgreen DVI VARI Carte2 Carte4 TVI VARI Vogelmann3 GNDVI Carte? G RVI 
0 SPVI2 Carte3 SPVI GRVI Carte3 VOG2 VOG2 Dattl VOG1 SPVI RVI NRI SR [7s2-690] 
OR LogR VR CR FD-0.3 FD-0.6 FD-0.9 FD-1.2 CWT-1 CWT-3 CWT-5 CWT-7 CWT-9 
光谱 变换 形式 
图 5 不 同 数学 变换 下 的 植被 指数 组 合 
Fig. 5 Vegetation index combination under different mathematical transformations 
表 2 棉花 叶片 叶绿素 含量 与 不 同 光谱 变换 下 多 植被 指数 的 相关 性 分 析 
Tab.2 Correlation analysis of chlorophyll content and vegetation index under different 
spectral transformations of cotton leaves 
光谱 变换 植被 指数 相关 性 范围 
OR SAVI,DVI,ND,RDVI,SPVI,SPVI2 0.24~0.26 
LogR VOG1,NDCI, Datt2,ND,RVI,NDVI,GNDVI, PSSRa, Clgreen, carte3 —0.27~-0.37 
FD-0.3 SAVI,OSARI, DVI, ND, RDVI, Carte3 , Carte4 , Carte2 ,SPVI,SPVI2 —0.28~0.29 
FD-0.6 DVI,SPVI2,SAVI,OSARI,EVI,ND,RDVI, Carte2, Carte4 , VOG2 —0.33~0.34 
FD-0.9 Vogelmann3 , RVI, DVI, SRyo7s-700;, Mndvins, ND, VOG1,NVI,TVI, VOG2 —0.23~0.30 
CWT-3 Clred, PSSRa, MSR», RVI, Lichtenthaler2 , VOG1 ,SIPI,Clgreen, GNDVI,SPVI —0.27~0.28 
CWT-7 REP, MCARLI, Vogelmann! , Datt!,GRVI,MTVI1 , cartel ,Ro,G,NRI -0.31~0.34 


Fy 0.858 和 0.894。 分 数 阶 变换 中 ,0.9 阶 的 RFR 模型 


注 :OR 为 原始 光谱 反射 率 ;LogR 为 对 数 光 谱 变换 ;FD-0.3 为 分 数 阶 0.3 阶 变换 ;CWT-3 为 小 波 变换 第 3 个 尺度 ,下 同 。 


量 的 重要 生化 参数 ,利用 光谱 技术 对 作物 的 叶绿素 


的 训练 集 尼 为 0.916 ,验证 集 尼 为 0.949 , 均 高 于 其 他 
阶 次 RFR 模型 的 尼 ,0.9 阶 的 WOA-RFR 模型 的 训练 
集 尼 为 0.920 ,验证 集 尼 为 0.955 , 比 原始 光谱 植被 指 
数组 合 WOA-RFR 模型 的 尼 提 高 了 4.5%。 在 小 波 
变换 尺度 上 ,第 7 尺度 上 的 WOA-RFR 模 型 的 训练 集 
尼 最 高 (CR=0.934) ,综合 分 析 , 上 述 模型 中 基于 分 数 
阶 0.9 阶 变换 上 的 多 植被 指数 组 合 建 模 模 型 WOA- 
RFR AY R E , RMSE 和 MRE 低 , 建 模 效 果 最 佳 ,可 以 
较为 准确 的 佑 测 棉 花 叶 片 叶绿素 含量 。 

从 预测 值 和 实测 值 的 拟 合 图 可 以 看 出 (图 6)， 
不 同 数学 变换 下 验证 集 的 模型 预测 效果 比 建 模 集 
的 预测 效果 好 ,基于 分 数 阶 0.9 阶 次 的 组 合 植被 指 
数 建立 的 模型 预测 效果 好 , 除 个 别 数据 点 其 他 都 较 
紧密 分 布 在 拟 合 线 周围 , 拟 合 效果 较 好 。 


3 讨论 


叶绿素 含量 一 直 是 评价 农作物 长 势 .产量 和 质 


含量 进行 定量 反 演 是 精准 农业 研究 的 热点 。 鉴 于 
此 ,本 文 以 南 疆 棉花 为 研究 对 象 , 基 于 不 同 的 光谱 
数据 变换 方法 计算 植被 指数 ,人 研究 结果 表明 ,数据 
变换 方法 能 够 提高 植被 指数 与 叶绿素 含量 的 相关 
性 ,这 与 前 人 的 研究 结果 一 致 , 童 新 等 ”基于 不 同 
阶 微分 高 光谱 植被 指数 的 牧区 草场 地 上 生物 量 
估算 研 究 中 表明 对 原始 高 光谱 反射 率 进行 微分 
处 理 , 有 助 于 光谱 数据 与 地 上 生物 量 的 相关 性 ; 
李 莉 婕 等 ”通过 高 光谱 定量 反 演 火龙 果 茎 枝叶 绿 
素 含 量 的 研究 得 出 数学 变换 和 连续 小 波 变 换 能 明 
显 提升 光谱 对 火龙 果 茎 枝叶 绿 素 含 量 的 估 测 能 

的 结论 ; 竟 霞 等 ”对 小 麦 条 锈病 遥感 监测 模型 构建 
的 研究 中 通过 对 原始 光谱 进行 分 数 阶 微分 处 理 , 提 
高 了 光谱 反射 率 与 小 麦 条 锈病 严重 度 的 相关 性 。 
另外 ,本 研究 结果 表明 ,利用 多 种 植被 指数 组 合 构建 
棉花 叶片 叶绿素 含量 的 模型 精度 较 其 他 研究 高 ,在 所 
有 光谱 数据 变换 上 RFR 和 WOA-RFR 模型 精度 都 高 
于 0.85, 而 在 分 数 阶 微分 的 0.9 阶 DVI ND, Vogel- 
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表 3 模型 建 模 结果 比较 
Tab. 3 Comparison of model modeling results 
光谱 变换 模型 算法 oe - = 

R RMSE MRE R RMSE MRE 

OR WOA-RFR 0.889 0.014 0.015 0.918 0.050 0.015 
RFR 0.860 1.307 0.019 0.862 1.240 0.018 

SVR 0.756 1.241 0.013 0.908 0.854 0.006 

LogR WOA-RFR 0.912 1.067 0.016 0.894 1.081 0.015 
RFR 0.884 1.367 0.020 0.858 1.280 0.020 

SVR 0.658 1.460 0.015 0.943 0.651 0.006 

FD-0.3 WOA-RFR 0.917 1.019 0.014 0.878 1.060 0.016 
RFR 0.872 1.357 0.020 0.850 1.269 0.019 

SVR 0.609 1.569 0.017 0.866 1.035 0.010 

FD-0.6 WOA-RFR 0.890 1.085 0.015 0.925 1.008 0.014 
RFR 0.872 1.379 0.019 0.846 1.388 0.021 

SVR 0.451 1.836 0.020 0.685 1.594 0.018 

FD-0.9 WOA-RFR 0.920 0.987 0.013 0.955 0.986 0.014 
RFR 0.916 1.250 0.018 0.949 1.207 0.017 

SVR 0.754 1.263 0.013 0.882 0.973 0.010 

CWT-3 WOA-RFR 0.895 1.091 0.015 0.922 1.086 0.014 
RFR 0.892 1.388 0.020 0.900 1.443 0.020 

SVR 0.779 1.172 0.010 0.906 0.835 0.007 

CWT-7 WOA-RFR 0.934 0.946 0.013 0.911 1.082 0.015 
RFR 0.913 1.286 0.019 0.841 1.281 0.020 

SVR 0.522 1.724 0.019 0.791 1.254 0.009 


预测 值 


预测 值 


注 : 尼 为 决定 系数 ;RMSE 为 均 方 根 误差 ;MRE 为 
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图 6 实测 值 与 预测 值 散 点 图 


Fig.6 Scatter plots of measured and predicted value 
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F 均 相对 误差 ;WOA 为 鲸鱼 优化 算法 ;RFR 为 随机 森林 回归 模型 ;SVR 为 支持 向 量 机 回归 模型 。 
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化 而 变化 ,因此 ,在 今后 研究 中 加 强 该 模型 在 不 同 的 长 绒 棉 叶绿素 含量 估算 中 ]. 激光 与 光电 子 学 进展 , 2022, 59 
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(2) 基于 不 同 光 谱 数 据 变换 的 不 同 植 被 指数 组 Sinan, Wang Quan, Jin Jia, et al. Application of hyperspectral indi- 
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Estimation of cotton leaf chlorophyll content based on combinations 
of multi-vegetation indices 


Areziguli ROZI™, Mamat SAWUT'**, HE Xugang™, YE Xiaowen™™ 
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Environment Modelling of Higher Education Institute, Urumqi 830017, Xinjiang, China) 


Abstract: Chlorophyll content is a crucial indicator for characterizing vegetation growth. In this study, we 
utilized high-spectral technology to rapidly monitor the chlorophyll contents of cotton leaves. We collected 125 
cotton leaf seedling samples from Xinjiang and measured their chlorophyll content and spectral data. To achieve 
this, we employed various spectral preprocessing techniques and used a combination of vegetation indices. 
Subsequently, we constructed a whale optimization algorithm/random forest regression (WOA-RFR) quantitative 
inversion model for cotton leaf chlorophyll content. Finally, we conducted a comparative analysis, contrasting the 
results of the WOA-RFR model with those obtained from the support vector regression (SVR) and RFR models. 
The results indicated that the spectral transformation methods (logarithm transformation, fractional order 
differentiation, and wavelet transformation) effectively improved the correlation between the vegetation indices 
and the chlorophyll content. We also found that the best inversion performance was achieved with the WOA-RFR 
model using a fractional order differentiation with a transformation order of 0.9 and the Vogelmann3, RVI, DVI, 
SR;s7s-720, Mndvins, ND, VOG1, NVI, TVI, VOG2 combined vegetation indices. The model exhibited R’ values of 
0.920 and 0.955 for the training set and validation set, respectively. The corresponding RMSE values were 0.987 
and 0.986, while the MRE values were 0.013 and 0.014. Compared to the RFR and SVR models, the WOA-RFR 
model demonstrated higher predictive accuracy, and the optimization effect of the WOA algorithm was evident. 
As a result, this study provides valuable decision- making support for accurately quantifying cotton leaf 
chlorophyll content. 


Keywords: combination of vegetation index; cotton; chlorophyll content; whale optimization algorithm 


